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High-precision calculations of hadron spectroscopy are a crucial task 
for Lattice QCD. State-of-the-art techniques are needed to disentangle the 
contributions from different energy states, such as solving the generalized 
eigenvalue problem (GEVP) for zero-momentum hadron correlators in an 
efficient way. We review the method and discuss its application in the 
determination of the _Bs-meson spectrum using (quenched) nonperturbative 
HQET at order l/mt. 

PACS numbers: 11.15.Ha 12.38. Gc 12.38. Qk 12.39.Hg 

1. Introduction 

Systematic errors in lattice simulations include not only finite-size and 
discretization effects, dependence on the chiral extrapolation and quenching, 
but also systematic effects due to contamination from higher excited states 
in the calculation of energy levels (masses) of quark bound states. This is 
because masses and decay constants are computed in lattice QCD from the 
exponential decay of Euclidean correlation functions C(t), which are built 
from composite fields with the quantum numbers of a given state (see e.g. 
[1]). More precisely, in the simulation, one evolves gluon fields (the link 
variables U) in the Monte Carlo dynamics associated with the partition 
function 



= J VUe-^^ deiK{U), (1) 
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where Sg is the gluonic action and K{U) is the Dirac operator. (The 
quenched approximation corresponds to det K = 1.) Although the quark 
fields ip are not evolved directly, information about them may be obtained 
once the link configuration is produced, since quark propagators are given 
by (i/jip) = (K^^) , i.e. by computing the inverse of K{U) for each configu- 
ration of the link variables and averaging over the configurations produced. 
Similarly, the computed inverse of K{U) may be used to build interpolators, 
i.e. products of creation and annihilation operators with the quantum num- 
bers of the desired bound state and a good overlap with the hadron wave 
function on the lattice. The result is then averaged over the configurations 
to yield the (Euclidean) correlators C{t) as 

C{t) = Y.{J{^,t)J{0,0)) = (O(t)O(O)), (2) 

X 

where J(x, t) = ipTip and T is the appropriate Dirac matrix (e.g. F = 75 
for pseudoscalar mesons). Finally, one determines masses and decay con- 
stants by identifying the various contributions to the spectral decomposition 

00 

(O(i)O(O)) = Cit) = Y: |(n|O|0)pe-^"S (3) 

n=l 

where \n) are eigenstates of the Hamiltonian (i.e. the logarithm of the trans- 
fer matrix) and all energies have the vacuum energy subtracted. Also, 
we assume Hermitean operators and a large enough time extent of the lat- 
tice to yield the simple exponential form above. At large t we expect to 
observe a plateau in the "effective mass" 

Ef{t) = log[C(t)/C(t + 1)] ^ ^1 + C)(e-(^2-^^)*) , (4) 

in such a way that the true ground-state energy Ei may be estimated. 

Clearly, determining masses and decay constants from the correlators 
in Eq. ([3]) is not an easy task. To see this, consider the first correction 
above, given by exp[— (£'2 — ^i)^] with a positive coefficient, in the case 
of a heavy-light system. For typical energy differences of a few hundred 
MeV, a plateau can only be reached for t around 1 fm, but by then the 
signal has started to compete with noise, even for improved static-quark 
discretizations. This is a general problem and trying to determine sub- 
leading corrections by multiple-exponential fits leads to large systematic 
errors. One alternative is to use more sophisticated fitting methods, such 
as Bayesian fitting, evolutionary fitting and NMR- inspired methods. An- 
other way to ensure better precision is inspired by the variational method 
in quantum mechanics and consists in increasing the basis of interpolators 
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to build a matrix of correlators Cij{t), for which a Generalized Eigenvalue 
Problem (GEVP) is formulated (see e.g. [2]). One then considers all-to-all 
propagators [3] instead of simple point sources as indicated in Eq. [2l 

The GEVP is a valuable tool to reduce systematic errors in the above 
determinations, and thus to deliver high-precision tests of QCD. We sum- 
marize the derivation of an optimal use of the method in Section [2] below 
and describe its application to spectrum calculations of Bg mesons in non- 
perturbative HQET to order l/m^ in Section O 



2. The Method 

The GEVP is defined by 

C{t)Vn{t,to) = Xn{t,to)C{to)Vn{t,to), (5) 

where t > to and C{t) is now a matrix of correlators, given by 

oo 

Cij{t) = {Oi{t)Oj{0)) = Y.e-^-'-^ni'fnj, i,j = l,...,N. (6) 

The chosen interpolators Oj are taken (hopefully) linearly independent, e.g. 
they may be built from smeared quark fields using N different smearing 
levels. The matrix elements are defined by 

'i'ni = (^n)i = {n\di\0) , (m|n) = 5mn ■ (7) 

One thus computes Cij for the interpolator basis Oi from the numerical 
simulation, then gets effective energy levels and estimates for the matrix 
elements ^'^j from the solution A„(i, io) of the GEVP at large t. For the 
energies 

it is shown 12J that Ef{t,to) converges exponentially as t ^ oo (and fixed 
to) to the true energy En- However, since the exponential falloff of higher 
contributions may be slow, it is also essential to study the convergence as 
a function of to in order to achieve the required efficiency for the method. 
This has been done in ^ , by explicit application of (ordinary) perturbation 
theory to a hypothetical truncated problem where only A'^ levels contribute. 
The solution in this case is exactly given by the true energies, and corrections 
due to the higher states are treated perturbatively. We get 

Ef{t,to) = En + en{t,to) (9) 
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Fig. 1. Schematic representation of energy levels, showing how a solution of the 
GEVP for conveniently chosen t, to yields conversion to the asymptotic state con- 
trolled by a much larger energy gap than usual (i.e. the one represented by the 
longer arrow instead of the shorter one) . 



for the energies and 



(10) 



for the eigenstates of the Hamiltonian, which may be estimated through 

Qf{t,to) = Rnid,Vnit,to)), (11) 



Rn = {Vnit,to), C{t)Vnit,to)) ^''^ 



An(to + g, tp) 
An (to + 2a, to) 



t/2 



(12) 



In our analysis we see that, due to cancellations of t- independent terms 
in the effective energy, the first-order corrections in en(t, to) are independent 
of io and very strongly suppressed at large t. We identify two regimes: 1) 
for to < t/2, the 2nd-order corrections dominate and their exponential sup- 
pression is given by the smallest energy gap \Em — En\ = AEm,n between 
level n and its neighboring levels m; and 2) for to > t/2, the Ist-order cor- 
rections dominate and the suppression is given by the large gap A£'iv+i,n- 
Amplitudes 7rn„/(t,to) get main contributions from the first-order correc- 
tions. For fixed t — to these are also suppressed with A£'7v+i,n- Clearly, the 
appearance of large energy gaps in the second regime improves convergence 
significantly. A pictorial illustration of the improvement in shown in Fig. [TJ 
We therefore work with t, to combinations in this regime. 

A very important step of our approach is to realize that the same pertur- 
bative analysis may be applied to get the leading corrections to correlators 
in an effective theory, such as corrections of order 1/m;, to the static case 
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in HQET correlation functions. These are given by 

C.,{t) = Cf\t) + uC]l'^'^{t) + 0{u') , (13) 

where the combined O^l/rrih) corrections are symboHzed by the expansion 
parameter w. FoUowing the same procedure as above, we get similar expo- 
nential suppressions (with the static energy gaps) for static and 0{l/mb) 
terms in the effective theory. We arrive at 

Ef{t,to) = Ef^''-\tM+^Ef'^''^^{tM) + 0{u^) (14) 

with 

'^*^*(t, to) = + PT' e-^^?^?i>" * + ..., (15) 

i?f'i/-Mt,io) = E}/^^ + [/^y-" - /3^*-4 A<;:\''^Je-^^-+V* + .(16) 

and similarly for matrix elements. 

An application of the methods described in this section is discussed next. 

3. Application to Nonperturbative HQET 

High-precision hadronic matrix elements are a key ingredient as theo- 
retical inputs in B physics and are ideally obtained from lattice-QCD simu- 
lations. However, currently used lattices are not large enough to represent 
simultaneously the low-energy scale of Aqcd , which requires a large physical 
lattice size, and the high-energy scale of the heavy-quark mass m^, which 
requires a very small lattice spacing a. A promising alternative is to con- 
sider (lattice) heavy-quark effective theory (HQET), which allows for an 
elegant theoretical treatment, with the possibility of fully nonperturbative 
renormalization [5|. 

HQET provides a valid low-momentum description for systems with one 
heavy quark, with manifest heavy-quark symmetry in the static limit. The 
heavy-quark flavor and spin symmetries are broken at finite values of rrib 
respectively by kinetic and spin terms, with first-order [i.e. 0{l/mi))] correc- 
tions to the static Lagrangian incorporated by an expansion of the statistical 
weight in 1/m;,, such that the symmetry-breaking operators are treated as 
insertions into static correlation functions. This guarantees the existence of 
a continuum limit, with results that are independent of the regularization, 
provided that the renormalization be done nonperturbatively. Masses and 
decay constants are expanded as sums of a static and an 0(l/mf,) contri- 
bution, in terms of the parameters of the effective theory and of the bare 
energies and matrix elements, which are computed in the numerical sim- 
ulation. The divergences (with inverse powers of a) in these parameters 
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are cancelled through the nonperturbative renormalization, which is based 
on matching the HQET parameters to QCD on lattices of small physical 
volume — where fine lattice spacings can be considered — and extrapolat- 
ing to a large volume by the step-scaling method. This analysis has been 
recently completed for the quenched case Using the computed HQET 
parameters, we have carried out a study [7] (see also [8]) of spectrum and 
decay constants for Bg mesons applying the GEVP method as described in 
the previous section. 

We have employed two lattice actions for the static quark, lattices of spa- 
tial extent L ~ 1.5 fm with three lattice spacings and all-to-all strange-quark 
propagators constructed from approximate low modes, with 100 configura- 
tions. The interpolating fields were obtained from a simple 7075 structure 
and 8 levels of Gaussian smearing for the strange-quark field. The resulting 
(8 X 8) correlation matrix may be conveniently truncated to an x one 
and the GEVP solved for each A^, so that results can be studied as a func- 
tion of N. We have used two methods for picking a basis from the above 
interpolators and checked that both yielded the same results. 

The combined use of nonperturbatively determined HQET parameters 
and efficient GEVP solution allowed us to reach a precision of a few percent 
in matrix elements and of a few MeV in energy levels, even with only a 
moderate number of configurations. A corresponding study for the Nf = 2 
case is in progress. 
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